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Abstract 

Density functional theory with local spin density approximation has been used to propose possible 
room temperature ferromagnetism in N-doped NaCl-type BaO. Pristine BaO is a wide bandgap 
semiconductor, however, N induces a large density of states at the Fermi level in the nonmagnetic 
state, which suggests magnetic instability within the Stoner mean field model. The spin-polarized 
calculations show that N-doped BaO is a true half- metal, where N has a large magnetic moment, 
which is mainly localized around the N atoms and a small polarization at the O sites is also 
observed. The origin of magnetism is linked to the electronic structure. The ferromagnetic(FM) 
and antiferromagnetic (AFM) coupling between the N atoms in BaO reveal that doping N atoms 
have a FM ground state, and the calculated transition temperature (Tc), within the Heisenberg 
mean field theory, theorizes possible room temperature FM in N-doped BaO. Nitrogen also induces 
ferromagnetism when doping occurs at surface O site and has a smaller defect formation energy 
than the bulk N-doped BaO. The magnetism of N-doped BaO is also compared with Co-doped 
BaO, and we believe that N has a greater potential for tuning magnetism in BaO than Co. 
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I. INTRODUCTION 


Inducing magnetism in non-magnetic oxides is one of the main areas of interest in the 
research community. Generally, many oxides are insulator and, without magnetic impurity 
atoms, non-magnetic. However, recent theoretical and experimental research of oxides has 
show ferromagnetism with defects or light elements, e.g., C, N, Li." In the past decade, 
density functional theory (DFT) has played a vital role in either proposing new magnetic 
materials or elucidating the origins of defects-driven magnetism in oxidesTo date, 
there are a few oxides that show room temperature (RT) ferromagnetism, e.g., ZnO, Sn02, 
111203, Ce02.— Light elements can not only be used to develop magnetism in oxides, but 
can also be used to stabilize intrinsic defects in the host materials.— Through the combined 
efforts of theoreticians and experimentalists, we are now certain that magnetism in oxides 
induced by non-magnetic impurities.-‘^‘ 22 ‘ 2 ^ Now it is known that magnetism develops in 
the nonmagnetic oxides where the impurity atoms have a hnite local magnetic moment and 
these local magnetic moments will interact with each other to form a net magnetic moment 
in the host material. Those materials in which magnetism is induced by doping nonmagnetic 
impurities, for which the substitutional atom can have a hnite magnetic moment and the 
2p-electrons of the doped atom plays an essential role in governing magnetism in the host 
material, are generally considered as magnetic materials. These nonmagnetic impurity 
atoms can also form impurity bands in the bandgap of host material, and ferromagnetic 
behavior can be expected if the Fermi energy lies within these impurity bands. 

BaO is an oxide with interesting structural and electronic properties. It can be used as a 
NO3 storage device for catalysis.— BaO is also considered as a precursor to the well known 
ferroelectric perovskite oxide BaTi03, which can have either Ti02 or BaO terminated sur¬ 
face when BaTi03 is grown on a suitable substrate as a thin him. Possible magnetism in 
Ti02.- has been extensively studied, however, less attention has been given to impurities 
in BaO and there are no detailed theoretical and experimental studies on inducing mag¬ 
netism in it. It has been experimentally observed that bulk BaO naturally occurs in a B 1 
(NaCl) structure.— Recently, experimental reports have claimed the growth of ultrathin BaO 
hlms on SrTiO3(001) substrate, and the formation of BaO nanoparticles on reconstructed 
SrTiO3(001), while a locally ordered c (4 x 4 ) BaO structure is observed on the disordered 
sample surface.— Tan et.al,— also observed RT ferromagnetism in N-doped BaTi03 and 
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the origin of the magnetism was correlated with the presence of N. Therefore, studying the 
electronic properties of BaO is of great importance for the development of a new applica¬ 
tions, specially in the area of magnetism. Hence, we propose to add a new functionality, 
i.e. magnetism, in BaO. We show that N in BaO has a ferromagnetic (FM) ground state 
and the transition temperature Tc, which is the temperature at which a material goes from 
a paramagnetic(disorder phase)state to a magnetic phase(order phase) state, is above room 
temperature. 

II. COMPUTATIONAL METHODS 

To study the magnetism of N-doped BaO, we performed calculations in the framework of 
density functional theory, — using linear combination of atomic orbital basis as implemented 
in the SIESTA code — . A double-^ polarized basis set for all atoms was used. The local 
spin density approximation: 22 (LSDA) was adopted for describing exchange-correlation inter¬ 
actions. We used standard norm-conserving pseudopotentials — in their fully nonlocal form 
— . A cutoff energy of 400 Ry for the real-space grid was adopted. The sampling of fc-space 
was performed with Monkhorst and Pack (MP) scheme with a regularly spaced mesh of 
18 X 18 X 18 . Convergence with respect to fc-point sampling and cutoff energy was carefully 
checked. 

To investigate the magnetism and electronic structures of N-doped BaO, we considered 
the 2x2x1 (Ba403Ni), 2x2x2 (BasOrNi), 2 x 2 x 3(Bai20iiNi) supercells of the primitive 
unit cell of NaCl-type BaO. Calculations were also carried out using a 3 x 3 x 3 (Ba27 026 Ni) 
supercell. In all these supercells, N was doped at O site due to smaller difference in their 
atomic sizes and electro-negativities. For the magnetic properties of BaO (OOl)surface, we 
considered the conventional unit cell of BaO, and used different n x m x z symmetric slabs 
of thickness dehned as a monolayer of BaO. We investigated n = l,m = 1 , z = 8 (Surfi) 
and n = 2 ,m = 1 , z = 8 (Surf2) surfaces of BaO and doped N at O on both sides of the 
BaO surface. We added a vacuum region of about n-u 10 A, so that the two surfaces do not 
interact with each other through the vacuum region. Additional simulations, using plane 
waves plus pseudopotentials as implemented in the Quantum Espresso (QE) code^, were 
also carried out to further test the validity of our results. Atomic positions were relaxed, 
using conjugate-gradient algorithm,— until the residual Hellmann-Feynman force on single 
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atom converges to less then 0.05 eV/A. Test calculations, using LSDA+t/, were also carried 
out by considering the on-site Coulomb correction {U = 6.0 eV, our previously optimized 
value^^) between the p-orbital electrons of O.— 


III. RESULTS AND DISCUSSIONS 

First, we calculated the optimized lattice constant of NaCl-type BaO, which was found 
to be 5.40A , as shown in Fig. [I^a). The calculated lattice constant is comparable to the 
previous calculated 5.47A,— and experimental 5.52A— values. Our DFT estimated value is 
smaller than the experimental value due to the underestimation problem of DFT-LDA. Using 
the optimized lattice constant, we calculated the band structure, which is shown in Fig. [Dj^b). 
We see that BaO is a wide bandgap semiconductor and the calculated bandgap at the X 
point is ~ 1.62 eV. The band is degenerate at G point [27r/a(0, 0, 0)], but this degeneracy is 
removed at X point where the band has less curvature (large effective mass) compared with 
point G. The electronic density of states (DOS) shows that bands near the Fermi energy 
(Up) are mainly derived from the O p-orbitals [see Fig. [2][a)]. The LSDA calculated indirect 
(F-A)bandgap is ~ 2.0 eV, in agreement with the previous LDA calculation.— However, the 
experimental and GW calculated value is 3.88 eV— . We then used LSDA-|-[/, and our 
LSDA-I-U calculated bandgap is ~ 3.6, which is comparable with the GW calculated and 
experimental values. 

It is important to consider atomic relaxation in N-doped BaO. Usually, the doping im¬ 
purity at the host site can either compress or elongate the bond length depending on the 
atomic size of the impurity atom. In some cases the impurity atom can also be located at an 
interstitial site and form a defect complex.— We therefore relaxed all the atomic coordinates 
of BaO:N and no signihcant changes in the atomic position of N at O sites were found. To 
further confirm this, both BagOyNi and Bai 20 iiNi were considered, and calculated their 
optimized lattice constants, which are 5.40A for both the systems [see Fig. [11(a)]. This 
lattice constant is approximately equal to our optimized lattice constant of pristine BaO, 
which shows that N can easily dope the O sites without any structural distortion. 

Figure [2](a) shows the calculated total and atom projected (P) DOS of N-doped BaO in 
the non-magnetic state. For comparison purposes, the total PDOS of pristine BaO is also 
shown. It can be clearly seen that the DOS of doped and pristine BaO are the same in 
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the conduction and valance bands. However, near the Fermi energy {Ep) the DOS of the 
doped system is drastically changed and the material shows metallic behavior, in contrast 
to pristine BaO, which is an insulator. The PDOS demonstrates that the main peak at the 
Fermi energy arises mainly from contributions from the p-orbitals of N. The N atom in BaO 
forms an impurity band in the bandgap of BaO. A small impurity-derived peak can also 
be seen at 2.0 eV just below Ep, and this peak hybridizes with the p-orbitals of O atoms. 
Therefore, the O and N atoms are only hybridize near the Fermi energy. No signihcant 
changes in the PDOS of O can be seen when the O sites are doped with N, except a small 
shoulder, derived by the N p-orbitals, near the Fermi energy. The calculated results show 
that the top of the valence band consists of N 2p electrons, which are higher than O 2p 
electrons. The large DOS at the Fermi energy D{Ep) in the N-doped BaO shows that there 
is an instability towards magnetism within the Stoner mean held theory of magnetism.— 
Within the Stoner model, which was mainly proposed for itinerant electron system, the 
large D{Ef) can lead to a large Pauli susceptibility which is large enough for the band to 
split spontaneously, and magnetism in N-doped BaO can be expected. 

To conhrm that N-doped BaO can present magnetism within the Stoner mean held the¬ 
ory, we further carried out spin-polarized calculations. The spin-polarized electronic density 
of state of Ba 403 Ni is shown in Figj2]j^b). The total DOS clearly shows that N induces 
magnetism in BaO when doped at O sites. The spin-polarized structure shows that the 
non-magnetic impurity band in the bandgap [see Figj2][a)] is mainly due to the contribution 
from minority spins. Large exchange splitting at N site can also be seen, as expected in the 
Stoner mean held theory. For comparison, we also show the total DOS of pristine BaO and 
one can easily judge that the conduction band of Ba 403 Ni in the spin-up state is identical 
to that of BaO. The 2p states of nitrogen are located at the top of the valence band, which 
are mainly derived from the anion p-states, resulting in the hybridization of electron wave 
functions of the 2p orbitals. The majority spin-states of nitrogen also hybridize with oxygen 
2p states, promoting the the impurity band which connects to the top of the valence band. 
Such chemical process rise to a bandgap in the majority spins states. The minority spins 
on the other hand create an impurity band in the bandgap which includes Ep. As a result 
this impurity band is broadened having about 0.38 eV half-width, and the approximated 
exchange splitting is 0.80 eV. The majority spins are completely occupied and behaving 
just like an insulator, whereas the minority spins are partially occupied and have a metallic 
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nature. Such electronic structure is the hngerprint of a half-metal ferromagnet,— which 
has applications in the area of spin electronics. The total and local magnetic moments of 
Ba 403 Ni were also calculated, and the calculated total magnetic moment is 1.0 /ib per unit 
cell. The PDOS clearly shows that the magnetic moment arises mainly from contributions 
from the N p-orbitals and the local magnetic moment at the N site is ~ 0.74 /xb- A small 
induced magnetic moment (~ 0.20 /xb) at the nearby O site is also observed. The PDOS of 
O shows some unoccupied minority spin states, which were occupied in the pristine BaO, 
suggesting that the holes induced by N also localized at the O p-orbitals. We repeated 
the same calculations by doping N at O site in 2 x 2 x 2 and 2x2x3 supercells, and 
interestingly we found half-metallic behavior in all these doped systems [see Fig. |2][d and 
e)] and the magnetic moments were mainly localized at the N sites. It is also noticeable 
that N does not polarize the whole valance band, it polarizes only the valance electronic 
states near the Fermi energy. Such behavior of N in BaO is different from that of tran¬ 
sition metals in oxides/semiconductors, which usually polarize the whole valance band.— 
The band-width of N decreases and the band is narrowed, as the concentration of N is 
reduced in BaO [see Fig. |2][e)], which further increases the localization of 2p states of N. 
This localization of the N 2p states may further increase the observation/suppression of 
ferromagnetism/antiferromagnetism in BaO, which is discussed in the following paragraph. 
The total magnetic moment per N atom is 1.0 /xb, which indicates that the total magnetic 
moment is independent of the N concentrations. In all these N-doped BaO systems, the 
calculated magnetic moment is consistent with Bund’s rules indicating that the N dopant 
exists as a anion illustrating that each N impurity introduces one 2p hole and 

the strong Hund-type exchange makes the nitrogen 2p spin-up orbitals completely occu¬ 
pied, and the spin-down states are partially hlled [Figj2][b-e)]. As a result, one can expect 
the (ferro)magnetism in N-doped BaO to be stabilized mainly by the predominant double¬ 
exchange mechanism.— Before we address the question of FM or antiferromagnetic (AFM) 
coupling between the N atoms, the above calculations were also performed for Ba 27 026Ni 
and Ba3203iNi, and we conhrmed that N in BaO has spin-polarized band structure. We 
also conhrmed the above conclusions using the QE code. It is also important to comment 
on the magnetism of N-doped BaO using LSDA-l-17. The LSDA-l-17 calculated total and 
PDOS of Ba 403 Ni are shown in Fig. [2](c). As expected, the N forms an impurity band in 
the bandgap of BaO and N-driven oxygen p states can also be seen in the gap region of BaO. 
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Similar observations have also been observed in self-interaction correction (SIC) for C-doped 
BaO.^^ The local magnetic moment of N (0)is increased (decreased), bnt the total magnetic 
moment of the nnit cell remains nnchanged, i.e., l.O/iB- Similar conclusions were also drawn 
when LSDA-|-f/ calculations were performed for BagOyNi, and Bai20iiNi consistent with 
our previous work.— 

As predicted, N can induce half-metallic magnetism in insulating BaO. However, there 
remains a question on whether N-doped BaO can have FM order above room temperature, 
and whether the exchange coupling sufficiently large to give rise to possible RT ferromag¬ 
netism. Magnetism is a cooperative phenomenon and a single N in BaO can not determine 
the true magnetic ground state. Therefore, we consider 2x2x3 and 3x3x3 super¬ 
cells of BaO and doped two N atoms at different O sites. The distance between the two N 
atoms d was varied, as it was the FM and AFM coupling between them. For comparison 
purposes, we kept the same N-N separation in both supercells, and the results are summa¬ 
rized in Table [H It is enchanting to see that both BaO:N systems have FM ground state, 
and the total magnetic moment of the unit cell is 2.0/ie (i.e., 1.0 /tb per N). All the sys¬ 
tems remained half-metallic as well. The strength of the exchange interactions J can be 
judged from AE = Eafm — Efm, where Eafm{Efm) is the total energy of the supercell 
in the AFM (FM) state. This exchange energy can be further used to estimate Tq- Using 
the Heisenberg mean held model {ksTc = 2 AU/3),— the estimated Tc is found to be 
close to room temperature, indicating the possibility of room temperature ferromagnetism 
in BaO:N systems. One can also see that both the systems have different transition tem¬ 
peratures due to different N concentrations in BaO. Note that, usually, mean held theory 
overestimates Tc but it can give an indication of possible RT FM. Therefore, RT ferromag¬ 
netism in N-doped BaO is expected, if properly prepared. It is encouraging to mention that 
our estimated Tc matches with N-doped BaTiO^^ and other magnetic oxides^*^ indicating 
that the magnetism is mainly attributed to N impurities. 

The formation enthalpy Hj of BaO is calculated using Hj = U(BaO) — [E(Ba) -|- |E(02)], 
where E(BaO), E(Ba), and E{02) are the total energies of NaCl-type BaO, BCC Ba, 
and oxygen molecule, respectively. The calculated Hj was found to be —6.49 eV, which 
comparable with the previous calculated (-5.64^, -5.19^), and experimental (-5.74^) values. 
We followed our previous approaches and calculated the defect formation energy Ej oi N 
in BaO under Ba-rich and 0-rich conditions because the defect formation energy strongly 


7 


depends on the growth condition and chemical potential.— The calculated Ef of Ba 403 Ni, 
BasOyNi, and Bai20iiNi in Ba-rich (0-rich) conditions are -1.15 (5.33), -1.14 (5.34), and 
-1.15 (5.45) eV, respectively. It is clear to see that N in BaO has a negative defect formation 
energy under Ba-rich condition, which further increases the possibility of introducing N in 
BaO. Note that LSDA-|-t/ calculated Ef of Ba 403 Ni is about 6.6 eV in 0-rich condition, 
which is increased by including U parameter-similar behvior has also been observed in other 
oxides.— We then varied the distance d between the two N atoms, using 2x2x3 and 3x3x3 
supercells, and calculated Ef under different growth conditions. The calculated Ef per N 
atom of both systems in the FM states is shown in Table [B The calculated Ef values are in 
the range of a single N doped BaO, and therefore we can not expect clustering of N in BaO 
when doping the O sites. The formation energy depends on the N concentration in BaO. We 
may safely conclude that all the studied systems have negative formation energies under Ba- 
rich environments. Such calculated results can help experimentalists to dope N in BaO. The 
role of intrinsic defects (vacancies) can not be ignored when a sample is synthesised, these 
intrinsic defects may interact with the impurity atom, and such interaction may further 
decrease Ef of the impurity atom.— There is also an experimental report on N-induced 
ferromagnetism in BaTi 03 ,— which suggests the possibility of N-doping at the O sites in 
BaO. Experimental work will be very helpful to further analyse the possibilities of N-doped 
BaO. 

To propose BaO for practical practical applications, e.g., in the area of thin Elms, it is 
essential to study the electronic and magnetic properties of N doped at surface O site in 
BaO, Surfi. We found that N can also induce magnetism in BaO (001) surface, and the 
magnetism is mainly localized at the N site as shown in Figl3](a,b). Surface oxygen atoms 
have negligible magnetic moment and N has a large magnetic moment (O.OT/ib) when doped 
at surface O sites, as expected for surface N atoms. No induced magnetic moment at the 
subsurface atoms were found, which is different from other doped oxides.->^ The PDOS of 
surface O atoms are shifted towards the Fermi energy compared to clean surface O atoms, 
indicating that the band gap of BaO(OOl) surface is reduced due to surface states driven by 
N-0 hybridization. The majority spins of N are occupied and the minority spins are partially 
occupied. The calculated surface spin densities of N, O, and Ba atoms are shown in Figl3](c), 
and it is evident to see no spin-polarization at the surface O and Ba atoms is observed. The 
spin density is mainly localized around the N atoms. We then used nx mx z surface (Surf 2 ) 
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and doped two N atoms at two different O sites and considered the FM and AFM conpling 
between the N atoms. Here, we have also fonnd that the FM conpling between the snrface 
N atoms has lower energy than that of AFM conpling. As the concentration of N atom is 
increases, a small magnetic moments at the snrface O atoms were also indnced, bnt the snb- 
snrface O atoms do not show any indnced magnetic moment [FigJHl^b)]. No snrface relaxation 
or reconstruction was observed when all the atoms were fully relaxed. The calculated Ef 
of N at surface O in Ba-rich (0-rich) condition is —1.58 (4.90) eV for Surfi. For Surf 2 , 
when the two N atoms were doped at the surface O sites, the Ef per N in Ba-rich (0-rich) 
conditions is —1.48 (5.00) eV when the distance between the two surface N atom is 3.82A, 
whereas when the separation between the two N atoms is 5.40A, the calculated Ef per N in 
Ba-rich (0-rich) condition is —1.60 (4.88) eV. Surface N-doped BaO has a smaller defect 
formation energy than the bulk N-doped BaO, similar to other oxides.— Note that a smaller 
surface Ef corresponds to a larger N-N bond distance suggesting that N prefers isolation 
than clustering in BaO(OOl) surface. 

We have shown that N in BaO can display possible room temperature ferromagnetism 
either in bulk or surface doped system, and hnally we would like to compare it with the 
Co-doped BaO. Our preliminary results clearly indicate that a Co atom at either O or Ba 
site induces magnetism with a larger(about 3 times) Ef than BaO:N systems. Therefore, 
we believe that N doping in BaO has the potential to easily tune magnetism and may be 
superior to Co doped BaO. Further experimental work would be needed to demonstrate the 
potential of BaO:N for possible applications in the area of spin electronics. We also believe 
that our prediction can also be helpful for BaO terminated surface of BaTiOa. 


IV. SUMMARY 

To summarize, DFT was used to investigate potential magnetism in N-doped NaCl-type 
BaO. Pristine BaO is known to be non-magnetic wide bandgap semiconductor, and N at 
O sites induce an impurity band in the bandgap of BaO. It has been discussed that N 
in BaO can also induce magnetism and the magnetism arises mainly from contributions 
by the N atoms. The minority impurity-driven band is narrowed with decreasing the N 
concentrations in BaO. The origin of magnetism was elucidated from the atom projected 
density of states. All the considered scenarios with different N concentrations in BaO showed 
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half-metallic behavior. The ferromagnetic and antiferomagnetic interactions between the N 
atoms were also considered, and they revealed that N-doped BaO is a ferromagnetic material. 
The transition temperatnre was also calcnlated using Heisenberg mean held model, and it 
was found to be above room temperature. The magnetic interaction between the N atoms 
is limited to the next-nearest-neighbor interactions. The surface magnetism of N-doped 
BaO(OOl) was also studied by doping surface O sites with N. N atoms at surface O sites 
also showed ferromagnetic ground state, and the magnetism was mainly limited to the BaO 
surface. No induced magnetic moments were observed at the sub-surface O sites. The 
defect formation energy showed that N has a smaller formation energy than the bulk N 
doped BaO. The magnetism of N-doped BaO was also compared to that Co-doped BaO, 
and it was concluded N-doped BaO can easily be achieved. LSDA-|-[/ calculations also 
supported ferromagnetism in N-doped BaO. 

The author acknowledges insightful discussions with Victor M. Garcia-Suarez and the 
cluster facilities of NCP, Pakistan. 
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Energy (eV) 


TABLE I. Calculated total magnetic moment M (in hb), AE = Eafm — EpM per cell (in eV), 
transition temperature Tq (in K), and defect formation energy Ej per N (in eV) of N-doped BaO 
(3 X 3 X 3). Babich and O^ich shows Ej under Ba-rich and 0-rich conditions, respectively. The 
first column shows N-N separation d (in A). Values in the brackets correspond to the 2x2x3 
supercell. 


d 

M 

AE 

Tc 

Ef 






Bdrich 

Or'ic/i 

3.82 

6.61 

2.0(2.0) 

2.0(2.0) 

0.064(0.077) 

0.035(0.089) 

495(597) 

271(689) 

-1.05(-1.07) 

-1.02(-1.08) 

5.42(5.40) 

5.45(5.40) 


(a) 



(b) 



L G X W K G 


FIG. 1. 


(a) The calculated lattice constant vs energy of pristine BaO (black line), BasOyNi (blue 


line), and Bai20iiNi( green line), (b) The band structure of pristine BaO, and the horizontal line 


(black) shows the Fermi energy, which is set to zero. 
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FIG. 2. The calculated total and atom projected density of states (PDOS) of Ba 403 Ni in the non¬ 
magnetic (a), where the green dotted lines show the PDOS(O-pris) and total DOS (pris) of pristine 
BaO. The spin-polarized (magnetic) total and PDOS of Ba 403 Ni, BasOrNi, and Bai20iiNi are 
shown in (b), (d), and (e), respectively. Panel (c) shows LSDA-|-[/ DOS of Ba 403 Ni. The green 
dotted lines in the total DOS panels show the total DOS of the pristine BaO (pris) of the same 
supercells. The total magnetic moments per N atom and local magnetic moments of N and O are 
also shown. The O and N PDOS are represented by solid (red) lines. The vertical lines represent 
the Fermi energy, which is set to zero. 
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FIG. 3. The LSDA calculated surface atomic projected DOS of Surfi(a) and Surf 2 (b). The top 
panels show the PDOS of surface O atoms (red lines) along with the PDOS of clean surface O 
atoms (green lines). The middle panels show the PDOS of surface N atoms (red lines). The bottom 
panels show the PDOS of sub-surface O atoms denoted as Osu- The local magnetic moments of 
N and O atoms are also shown. The vertical lines show the Fermi energy, which is set to zero. 
The calculated surface spin density of Surfi(c), where red, green, and yellow balls represent Ba, O, 
and N atoms, respectively. The spin density of Surfi is calculated on a large surface to show zero 
polarization at surface O atoms. The inset shows the scale of the spin density. 
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